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Abstract I briefly describe a simple routine for emission-) 
line profiles fitting by Gaussian curves or Gauss- 
Hcrmite series. The PROFIT (line-PROfilc FlTting) rou- 
tine represent a new alternative for use in fits data 
cubes, as the ones from Integral Field Spectroscopy or 
Fabry-Perot Interferometry, and may be useful to bet- 
ter study the emission-line flux distributions and gas 
kinematics in distinct astrophysical objects, such as 
the central regions of galaxies and star forming regions. 
The profit routine is written in IDL language and is 
available at |http: / / www.ufsm.br / rogemar / software. html | 

The profit routine was used to fit the [Fen] A = 
1.257/im emission-line profiles for about 1800 spectra 
of the inner 350 pc of the Seyfert galaxy Mrk 1066 
obtained with Gemini NIFS and shows that the line 
profiles are better reproduced by Gauss-Hermite se- 
ries than by the commonly used Gaussian curves. The 
two-dimensional map of the Gauss-Hermite moment 
shows its highest absolute values in regions close to the 
edge of the radio structure. These high values may 
be originated in an biconical outflowing gas associated 
with the radio jet - previously observed in the optical 
[O in] emission. The analysis of this kinematic com- 
ponent indicates that the radio jet leaves the center of 
the galaxy with the north-west side slightly oriented to- 
wards us and the south-east side away from us, being 
partially hidden by the disc of the galaxy. 
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1 Introduction 

Integral Field Spectroscopy (IFS) and Fabry-Perot In- 
terferometry are powerful tools to do a two-dimensional 
analysis of the physical properties of several types of 
astronomi cal objects, such as the central region of nor- 
mal (e.g. iRodrigues et al.l [20091; IPeletier et al.l 12007 : 



ies (e.g. iDavies et al.l 120091 120071 Ifficks et al.l 12009 : 
Riffel Stqrchi-Bergmann fc Naearl l20ld IRiffel et al 



2010 



Emsellem et al. l2007t iDiaz et a I2OO6I) and active gal 



lax- 



lallbl. 120081 l2006t IStorchi-Bergmann et al.l 12009 . 
2007t iFathi et al l 120061). star forming re gions 



(e.g. lBlum fc McGregorll2009l:lBarbosa et al.ll2008h and 
young stellar objects (e.g. [Beck et al .l2008llMcGregor et al| 



120071: iTakami et al J 120071 ). The final result of the data 
reduction of the above techniques is a data cube con- 
taining hundreds to thousands individual spectra to 
be analyzed, thus a common problem among the stud- 
ies cited above is how to extract the information from 
these data cubes. A manually inspection of each spec- 
trum is an exhaustive task and demand much time, 
so automated methods are needed to properly measure 
physical parameters from the data cubes. 

The most common method used to study the gaseous 
distribution and kinematics is based on the fitting of 
the emission-line profiles by Gaussian curves. Nev- 
ertheless, it is commonly reported in the literature 
the presence of asymmetries in the emission-line pro- 
files, which are not well represented by Gaussian 
curves (e.g 



Barbosa et al. 2009; Komossa et al.l 2008 



Riffel et al 



2009at IRiffel fc Storchi-Bergmannl boiol) 



Sev eral methods have been developed to fit line profiles 
(e.g ISarzi et al 20061 and IRAF fitprofs and SPLOT 
tasks), but most of these methods are based on the 
fitting of Gaussian (or Lorentzian) functions, in which 
the information of the wings of the line profile can be 
lost. A recent developed method to extract information 
from data cubes is the PCA tomography, which uses 
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Principal Component Analysis (PCA) to transform the 
system of correlated coordinates into a system of un- 
corrected coordinates o rdered by principal components 
of decreasing variance ( Steiner et al. 20091) . Neverthe- 
less, in some cases the 'traditional' line-profile fitting 
method must be additionally used to properly extract 
the flux distribution and kinematics of the emitting gas. 

In this work I present an automated line-PROfile 
FiTting routine (profit) to be used to extract the 
gaseous kinematics and flux distribution from fits data 
cubes. This routine is written in IDI0 language and 
allows the fit of the observed profiles by Gauss-Hermite 
series or Gaussian curves. The fit of Gauss-Hermite 
series has been chosen because it preserve the veloc- 
ity information of the emitting gas by the fitting of 
the wing of the emission-line profiles. Such informa- 
tion could be lost in the fit of a single Gaussian curve 
for an asymmetric emission-line profile. Another ad- 
vantage of the Gauss-Hermite profile is that it can be 
easily be implemented in an automated routine than 
multiple Gaussian fit - which would also preserve the 
velocity information. 

The paper is organized as follows. In Section 2 I 
present the formalism of the Gauss-Hermite and Gaus- 
sian functions; Sec. 3 presents the profit routine and 
in Sec. 4 I discuss an application of the routine for the 
Seyfcrt galaxy Mrkl066. Sec. 5 presents the final re- 
marks of the present work. 



2 Gauss-Hermite versus Gaussian 



The Gauss-Hermite series can be written as (e.g. 




Fig. 1 Comparisson of Gaussian curves (dotted lines) with 
Gauss-Hermite series (continuous lines) for the hz and /14 
values shown at the top-left corner of each panel. The am- 
plitude, central wavelength and a are the same for both 
functions and have arbitrary values. 



If the emission-line profile is similar to a Gaus- 
sian we can truncate the sum on n = 4 and assume 
hp = Hq(w) = 1, hi = hi = (jvan der Marel fe Franx 
19931 ). This is a good approximation if the emission- 



van der Marel fe Franxll993tlGerhardl99dlCappellari fe 
2004): 



line presents an asymmetric profile, such as the blue 
or red wings frequently observed for the line emis- 
sion from ion ized gas in the narrow line region of ac- 
tive galaxies Jpaffel. Storchi-Bergmann fe Nagar 2010 ; 
IMM(atijH l2009al iKomossa et al.ll2008h . Usingthe ap- 
proximation above, the Eq.fTJcan be written as 



GH = 



Aa(w) 



where 

A -A, 



and 

a(w) 



>2tt 



(1) 



(2) 



(3) 



A is the amplitude of the Gauss-Hermite series, A c is the 
peak wavelength, hj are the Gauss-Hermite moments 
and Hj(w) are the Hermitc polynomials. 

1 http:/ /www. ittvis.com/ 



GH 
where 
H 3 (w) 

and 



Aa(w) 



[1 + h 3 H 3 (w) + h 4 H 4 (w)} , 



V6 



(4) 



(5) 



(4w 4 - 12i 



3). 



(6) 



The /13 Gauss-Hermite moment measures asymmet- 
ric deviation from a Gaussian profile, such as blue or 
red wings, while the hi moment quantify the peaki- 
ness of the profile, with h^ > for a more peaked 
and /14 < for a broader profile than a Gaussian 
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curve. A particular case of Eq. [4] is /13 = hi = 0, 
when it becomes a Gaussian curve. In Figure [TJ I 
present a sample of profiles for Gauss-Hcrmite series 
with distinct /13 and /14 values (continuous lines) and 
Gaussian curves (dotted lines). The amplitude, cen- 
tral wavelength and a are the same for Gaussian and 
Gauss-Hcrmite functions and have arbitrary values. 
The /13 and /14 moments are shown at the top-left 
corner of each panel and have values typically ob- 
serve d for the narrow line region of active galaxies 
(e.g. Riffel fc Storchi-Bergmann 1201(1 The compar- 
isson of observed e mission-line profiles, such as those 
of Mr k 1066 from iRiffel. Storchi-Bergmann h, Nagar 
( 2010l ) (see Sec.|4]), with the profiles shown in Figured] 
suggests that the observations are better reproduced by 
Gauss-Hermite series than by Gaussian curves in most 
cases. 



3 The profit routine 

The model of each emission-line profile is constructed 
by the sum of Eq. f?] with a linear equation, in or- 
der to represent the underlying continuum emission. 
The resulting equation contains seven free parameters 
(^4, A c , a, /13 /14 plus two parameters for the linear equa- 
tion), which can be determined by fitting the line pro- 
files. In case of Gaussian fitting the /13 and /14 are fixed 
at zero and thus the remaining 5 parameters may be 
obtained from the fit of the observed profile. These 
parameters can be obtained by solving a Least-squares 
problem. 

The profit routine was written in IDL language and 
performs the fit of the observed profile using the MPFI 
routine, which the is the M I NPAC K1 implementation 
( More. Garbow fc Hillstrom 1980l) of the Levenberg- 
Marquardt method for nonlinear least-squares prob- 
lems. The IDL language was chosen because it is exten- 
sively used in astronomy and allows read the data cube 
from standard fits format to an array using the NASA- 
Goddard Space Flight Center IDL Astronomy User's 
LibrarjU. The algorithm recovers the emission-line flux 
distribution and kinematics as follows: 

1- Initial guesses for the centroid wavelength and veloc- 
ity dispersion are given by the user. The initial guesses 
for /13 and /14 are fixed at zero; 

2- The input data cube in standard fits format (in which 
the spatial dimensions are in the x and y-axis and spec- 
tral pixels are in the z-axis) is converted into an array. 



2 The MPFIT routine can be obtained from the Markwardt IDL 
Library at http://cow.physics.wisc.edu/~craigm/idl/idl.html 

3 Thc IDL Astronomy User's Library can be obtained from 
http://idlastro.gsfc.nasa.gov 



Next steps are done individually for each spectrum: 

3- Calculates the spectral region to be fitted using the 
initial guess for the centroid wavelength and the spec- 
tral information (spectral sampling and initial wave- 
length) contained in the header of the data cube fits 
file. 

4- Normalizes the spectrum by its maximum value and 
obtain initial guesses for the parameters of the linear 
equation and amplitude of the Gauss-Hcrmite series (or 
Gaussian); 

5- Performs the nonlinear least-squares fitting of the 
observed profile by the adopted the model using the 
Levenberg-Marquardt method; 

6- Writes the solution to the output file; 

7- If \ 2 is l ess than a maximum value (defined by the 
user) the fitted parameters are used as initial guesses 
for the fit of the next spectrum. Otherwise uses the 
initial guesses of item 1; 

8- Repeats items 3-8 for all spectra. 

9- Writes the solutions to a Multiple Extensions FITS 
(MEF) file. The output MEF file will contain 7 exten- 
sions containing the: [0] emission-line flux distribution; 
[1] centroid velocity field; [2] velocity dispersion map; 
[3] /13 map; [4] h 4 map; [5] reduced x 2 map defined as 



E 



(O p - M v f 



1 



(N - A par ) ' 



(7) 



where O p is the observed spectra, M p is the best fit 
model, <Tq is the variance of the observed spectra, N 
is the number of spectral pixels (p) used in the fit and 
,/Vpar is the number of free parameters; and [6] flux dis- 
tribution obtained directly by integration the emission 
line profile and subtracting a continuum obtained by 
the average of continuum regions at both sides of the 
line profile. 



4 A First Application 

In Fig. [2] I present an example of the use of profit to 
fit the emission-line profile of [Fen] at A = 1.257 /im 
at 1" north-west of the nucleus of the Seyfert galaxy 
Mrk 1066. This spectrum was extracted from Gem- 
ini's Near-infrared Integral Field Spectrograph (NIFS) 

observations (program ID: GN -2008A-Q-30) within an 

apertu re of 0.3 x 0.3 arcsec 2 (see Riffel. Storchi-Bergmann fc Nag arj 
(2010) for a description of the observations and data 
reduction procedures). The observed profile is shown 
as a continuous line, the fitting of Gauss-Hcrmite series 
as a dashed line, the fitting of a single Gaussian as a 
dot-dashed line and the two Gaussian curves fit as a 
dotted line. The best Gauss-Hermite fit is obtained 



1.270 



1.272 



1 .274 



1.276 



Fig. 2 Example of fitting using profit routine for the 
[Fell]Af.2570 fim emission-line profile for a spectrum ex- 
tracted within an aperture of 0.3x0.3 arcsec 2 at 1 arcsec 
north-west of the nucleus of Mrk 1066. The continuous line 
represents the observed profile, the dashed line the result- 
ing Gauss-Hermite fit, the dot-dashed line the fit of a single 
Gaussian and the dotted line the resulting fit by two Gaus- 
sian curves. 



for a central wavelength A c = 12722.6 ± 0.2 A, a ve- 
locity dispersion of a = 114.3 ± 8.8kms~ 2 and higher 
order Gauss-Hermite moments /13 = —0.085 ± 0.017 
and h 4 = 0.145 ± 0.014 and resulted in a \ 2 = 0.014. 
The best fit obtained using a Gaussian curve has 
X 2 = 0.049 is obtained for A c = 12722.8 ± 0.2 1 and 
a = 106. 2± 11.9 kms~ 2 . Although the values obtained 
for A c and a from Gauss-Hermite and Gaussian func- 
tions are similar, Fig. [2] clearly shows that the observed 
profile is better reproduced by Gauss-Hermite series 
than by a Gaussian curve - the blue wing present on 
the observed profile is not reproduced by the Gaussian 
curve. This conclusion is also evidenced by the lower 
X 2 value obtained for the Gauss-Hermite fitting. In 
the other hand, the resulting \ 2 obtained for the fit- 
ting of the [Fe 11] emission- line profile by two Gaussian 
curves x\g ~ 0-011 is similar to those obtained from 
the Gauss-Hermite fitting. However, in an automated 
fitting routine it is hard to decide when an emission- 
line profile should be fitted by a single Gaussian and 
whe n it should be fitted by multip le Gaussian curves 
(e.g. Storchi-Bergmann et al. 2010[) . while for Gauss- 
Hermite series this decision is done simple by varying 

the h$ and hj moments. 

In Riffel fe Storchi-Bergmann d20ld ) we used profit 
to study the gaseous kinematics of the inner 700 x 700 pc 2 | 
of Mrk 1066 using Gemini NIFS observations. The pro- 
files of [P11] A1.1886/mi, [Fen] Al. 2570 fim, Pa/3 and 
H 2 A2.1218yum emission lines were fitted by Gauss- 
Hermite series in order to obtain the centroid velocity 



Fig. 3 Comparisson of the resulting \ 2 obtained from 
the fitting of Gauss-Hermite series (y-axis) and Gaussian 
curves (x-axis) to the [Fell] A1.2570/im emission-line profile 
of ~1800 spectra of the inner 700x700pc 2 of Mrk 1066. 



field, velocity dispersion map and /13 and /14 maps for 
ea ch emission line. The same obser vations presented 
in iRiffel fe Storchi-Bergmann] ( 2010 ) were used to fit 
the [Fell] Al. 2570 /xm by Gaussian curves in order to 
compare the resulting fit with the ones obtained using 
Gauss-Hermite scries. Figure [3] presents the compar- 
isson of the x 2 values from the Gauss-Hermite fitting 
(Xgh ~ y-axis) with the ones obtained for the fitting 
of Gaussian curves (xq - x-axis) for ^1800 spectra 
of the inner 700x700pc 2 of Mrk 1066 extracted within 
apertures of 0.05x0.05 arcsec 2 . As observed in this fig- 
ure, Xgh ^ s smaller than x% ^ or m °st spectra, indicat- 
ing that the [Fen] line profile in the central region of 
Mrk 1066 is better reproduced by Gauss-Hermite series 
than by Gaussian curves. A similar behavior is ob- 
served for [Pn], H2 and H emission lines at the same 
spatial region. 

In order to illustrate the importance of properly map 
the emission line profile wings Fig. [4] presents a map for 
the Gauss-Hermite moment obtained for the Pa/3 
emission in the central region of Mrk 1066. This map is 
simmilar to t he ones obtained for the [Fe n] an d [P 11] 
emission lines ( Riffel fe Storchi-Bergmannll2010l ) , show- 
ing several regions with values different than zero, indi- 
cating that the Pa/3 emission-line profile presents asym- 
metric deviation from a Gaussian curve, such as blue 
(negative values) and red (positive values) wings. In 
case of Gaussian fitting this information could be lost! 
The smallest values of up to —0.3 are observed at ~ 1" 
north-west of the nucleus in a regions close to the edge 
of the radio structure (thick black contours). Some 
high values are also observed near to the edge of the 
radio jet to south-east of the nucleus. The presence of 
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Fig. 4 Top: /13 Gauss-Hermite moment map for Pa/3 emis- 
sion line on the inner 700x700pc 2 of Mrkl066. The thick 
black contours are from the radio continuum emission of 
iNagar et all (|l999l ). Bottom: Ha Gauss-Hermite moment 
map. 



wings have also been observed bv lKnop et alj (|2001l ) for 
the near-IR emission lines using long slit spectroscopy 
along the PA=135°. These wings may be originated in 
an outflowing gas component driven by the radio jet, 
with the north-west side slightly oriented towards us 
and the south-east side away from us, being partially 
hidden by the disc of the galaxy. This interpretation is 
supported by the near-IR emission line kinematic maps, 
which show that the near-IR emitting gas presents at 
least three kinematics components: a rotating disk, an 
inflowing gas component and an o utflowing component 
dRiffel fe Storchi-Bergmai"rn1l2010h . The above interpre- 
tation is also in good agreement with optical observa- 
tions of the [O ill] emission, which seems to being orig- 
inated in a bi-cone oriented along the same position 
angle PAsal35° — approximately t he same orientation 
of the radio jet (jBower et al.lll995l ). 

The hi map, shown in the bottom panel of Fig. [4j 
presents values near to zero in most locations of the 
central region of Mrk 1066. However, some positive 
values are observed co-spatially with regions of lower ve- 
locity dispersion, indicating that part of the Pa/3 emis- 
sion originates from a colder gas than those which pro- 
duces most of the line emission. For more details on 
the gaseous kinematics of the inner 350 pc ra dius of 
Mrk 1066 see lRiffel fc Storchi-Bergmannl (|2010l ). 



5 Final Remarks 

I presented a new routine to fit emission-line pro- 
files from fits data cubes using Gauss-Hermite series 
or Gaussian curves, which provides a new alterna- 
tive to the study of the emission-line flux distribu- 
tion and kinematics for several astronomical objects. 
The profit routine is written in IDL and is available 
at |http: / / www.ufsm.br/rogemar/software.htmll The 
main advantages of profit compared with previous 
methods are: 

• It allows the fitting of the line profiles by Gauss- 
Hermite series as an alternative to the 'traditional' 
Gaussian curves used in most studies and thus bet- 
ter describe the kinematics of the emitting gas. 

• It is automated and can be applied directly on the 
final data cubes fits files from observations using most 
integral field units and Fabry-Perot interferometers. 

A first application has been discussed for the case 
of near-IR emission-line profiles from the inner 350 pc 
of the Seyfert galaxy Mrk 1066. The main scientific 
conclusions are: 



• The near-IR emission line profiles for this galaxy are 
better reproduced by Gauss-Hermite series than by 
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Gaussian curves, as indicated by the smaller % 2 ob- 
tained for the former. 
• The two-dimensional map for the /13 Gauss-Hermitc 
shows that the near-IR emission lines present asym- 
metric profiles, which can be explained as being orig- 
inated in an outflowing gas driven by the radio jet. 
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